x0 = [0; 0];
f0 = rosenbrock(x0(1), x0(2));
df0 = 
